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An  Extension  of  the  Argument  Principle  and  Nyquist 
Criterion  to  Systems  with  Unbounded  Generators 

Makan  Fardad  and  Bassam  Bamieh 


Abstract 

The  Nyquist  Stability  Criterion  is  generalized  to  systems  where  the  (open-loop)  system  has  infinite-dimensional 
input/output  spaces  and  a  (possibly)  unbounded  infinitesimal  generator.  This  is  done  through  use  of  the  perturbation 
determinant  and  an  extension  of  the  Argument  Principle  to  infinitesimal  generators  with  trace-class  resolvent. 


I.  Introduction 


The  Nyquist  criterion  is  of  particular  interest  in  system  analysis  as  it  offers  a  simple  visual  test  to  determine 
the  stability  of  a  closed-loop  system  for  a  family  of  feedback  gains  [1]  [2].  Extensions  of  the  Nyquist  stability 
criterion  exist  for  certain  classes  of  distributed  [3]  and  time  periodic  [4]  systems.  [3]  considers  distributed  systems 
in  which  the  open-loop  G{s)  belongs  to  the  algebra  of  matrix-valued  meromorphic  functions  of  finite  Euclidean 
dimension,  and  the  Nyquist  analysis  is  carried  out  by  performing  a  coprime  factorization  on  G(s). 

To  motivate  the  discussion  in  this  paper,  let  us  first  consider  a  finite-dimensional  (multi-input  multi-output)  LTI 
system  G(s)  placed  in  feedback  with  a  constant  gain  7/.  In  analyzing  the  closed-loop  stability  of  such  a  system, 
we  are  concerned  with  the  eigenvalues  in  C"*"  of  the  closed-loop  zl-matrix  If  s  is  an  eigenvalue  of  then 

it  satisfies  det[s/—  =  0.  Now  to  check  whether  the  equation  det[s/—  =  0  has  solutions  inside  C'*',  one 

can  apply  the  argument  principle  to  det[/  +  7G(s)]  as  s  traverses  some  path  £>  enclosing  C'*'.  More  concretely, 
let  us  assume  that  we  are  given  a  state-space  realization  of  the  open-loop  system.  Then  using 


det[/  +  7G(s)] 


det  [si  — 
det  [si  —  A] 


(1) 


if  one  knows  the  number  of  unstable  open-loop  poles  one  can  determine  the  number  of  unstable  closed-loop  poles 
by  looking  at  the  plot  of  det[/  +  7G(s)]|^gjj. 

But  in  the  case  of  distributed  systems  the  open-  and  closed-loop  infinitesimal  generators,  A  and  may  be 
operators  on  an  infinite-dimensional  Hilbert  space  X  and  can  be  unbounded.  Hence  it  is  not  clear  how  to  define  the 
characteristic  functions  det\sX  —  A\  and  det[sX  —  In  this  paper  we  find  an  analog  of  equation  (1)  applicable 
to  unbounded  A  and  A‘^^  and  use  operator  theoretic  arguments  to  relate  the  plot  of  det[X  +  to  the 

unstable  modes  of  the  open-loop  and  closed-loop  systems. 

Now  if  the  multiplicity  of  each  of  the  eigenvalues  of  A  is  finite  it  can  be  shown  that  det[X  +  7(/(s)]  is  still 
a  meromorphic  function  of  s  on  C,  and  one  may  be  tempted  to  use  the  methods  of  [3]  to  analyze  closed-loop 
stability.  But  if  the  open-loop  system  has  distributed  input  and  output  spaces,  then  [3]  requires  the  coprime 
factorization  of  an  infinite-dimensional  operator.  In  addition,  one  often  deals  with  systems  of  Partial  Differential 
Equations  (PDEs)  in  which  the  state-space  representation  is  the  natural  representation  and  it  is  more  convenient 
to  deal  directly  with  the  operators  A  and  rather  than  Q{s)  [see  example  in  Section  V]. 

Our  presentation  is  organized  as  follows:  We  lay  out  the  problem  setup  in  Section  II  and  describe  the  general 
conditions  for  stability  of  distributed  systems  in  Section  III.  Section  IV  contains  the  main  contributions  of  the 
paper  in  which  the  Argument  Principle  and  the  Nyquist  Stability  Criterion  are  extended  to  a  class  of  distributed 
systems.  The  theory  is  applied  to  a  simple  example  in  Section  V.  Proofs  and  technical  details  have  been  placed 
in  the  Appendix. 
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Notation 

S(r)  is  the  spectrum  of  T,  and  p{T)  its  resolvent  set.  an{T)  is  the  n*  singular-value  of  T.  ^{X)  denotes  the 
bounded  operators  on  the  Hilbert  space  X,  ^ao{X)  the  compact  operators  on  X,  and  ^i(X)  the  nuclear  (trace- 
class)  operators  on  X,  i.e.  operators  T  that  have  the  property  ^n(T)  <  oo;  ^i(X)  C  t^ooiX)  C  ^(X).  tr[T] 

denotes  the  trace  of  T  and  det[T]  its  determinant.  C’*'  and  C“  denote  the  closed  right-half  and  the  open  left-half 
of  the  complex  plane,  respectively,  and  j  :=  C(zo;^)  is  the  number  of  counter-clockwise  encirclements 

of  the  point  zq  G  C  by  the  closed  path  ip. 


II.  Problem  Setup 

Consider  the  open-loop  system  S° 

[dt^]{t)  =  [Ai^]{t)  +  [Bu]{t), 

y{t)  =  [CV'IW,  (2) 

where  t  G  [0,  oo)  with  the  following  assumptions.  The  (possibly  unbounded)  operator  A  is  defined  on  a  dense 
domain  ^  of  the  Hilbert  space  X  and  is  closed.  B  and  C  are  bounded  operators  on  X.  At  any  given  point  t  in 
time,  u,  y  and  i/;  belong  to  the  space  X  and  are  the  distributed  input,  output,  and  state  of  the  system,  respectively. 
We  will  refer  to  A  as  the  infinitesimal  generator  of  the  system.  We  may  also  refer  to  A,  B,  and  C  as  the  system 
operators.  The  open-loop  system  Sg  has  temporal  impulse  response  Q{t)  :=  C  B,  and  transfer  function 

g(s)  :=  C{sI-A)-^B.  (3) 

Next  we  place  the  system  S°  in  feedback  with  a  bounded  operator  yT,  |jX'||  =  1,  7  G  C.  This  forms  the 
closed-loop  system  shown  in  Figure  1  (left)  with  infinitesimal  generator  A'^^  :=  A  —  B  yT C.  We  separate  the 
function  T  from  the  gain  7  as  in  Figure  1  (right)  to  form  the  closed-loop  system  8“^*,  and  it  is  our  aim  here  to 
determine  the  stability  of  8“^*  as  the  feedback  gain  7  varies  in  C. 


Fig.  1.  Left:  The  spatially  periodic  closed-loop  system  as  the  feedback  interconnection  of  a  spatially  invariant  system  G  and  a  spatially 
periodic  multiplication  operator  F.  Right:  The  closed-loop  system  8“^*  in  the  standard  form  for  Nyquist  stability  analysis. 


We  also  make  the  following  assumptions: 

Assumption  (*):  There  exists  at  least  one  point  s  G  p{A)  such  that  {sX  —  A)~^  G  i^i(X), 

Assumption  (**):  p{A)  contains  a  right  sector  of  the  complex  plane  |  arg(z  —  a)  |  <  |  -f  (p,  (/?  >  0,  a  G  M. 

III.  Stability  of  Distributed  Linear  Systems 

A  semigroup  on  a  Hilbert  space  is  called  exponentially  stable  if  there  exist  constants  M  >  1  and  >  0 
such  that  ||e'^‘||  <  for  f  >  0.  It  is  well-known  [5]  [6]  that  if  A  is  an  ifinite-dimensional  operator,  then 

in  general  S(^)  C  C“  is  not  sufficient  to  guarantee  the  exponential  decay  of  ||e'^*||.  In  this  paper  we  focus  on 
systems  which  do  satisfy  the  so-called  spectrum-determined  growth  condition,  i.e.,  systems  for  which  S(,/l)  C  C“ 
does  imply  exponential  decay  of  the  semigroup.  Examples  of  such  semigroups  are  numerous  and  include  analytic 
semigroups  [7]  [8]. 

Thus  to  guarantee  the  exponential  stability  of  8‘^*,  it  is  necessary  and  sufficient  to  show  that  A^^^  =  A  —  B  yT C 
has  spectrum  only  inside  C“.  In  the  next  section  we  aim  to  develop  a  graphical  method  of  checking  whether 
or  not  C  C“.  Also,  henceforth  in  this  paper  wherever  we  use  the  term  stability  we  mean  exponential 

stability. 
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Fig.  2.  The  closed  contour  D  traversed  in  the  clockwise  direction  taken  as  the  Nyquist  path  as  r  — >  oo.  The  indentations  are  arbitrarily 
made  to  avoid  the  eigenvalues  of  A  (i.e.,  open-loop  modes)  on  the  imaginary  axis. 


IV.  The  Nyquist  Stability  Criterion  for  Distibuted  Systems 
A.  The  Determinant  Method 

As  discussed  in  the  Introduction,  we  aim  to  use  operator  theoretic  arguments  to  relate  the  plot  of  det[X  + 
to  the  unstable  modes  of  the  open-loop  and  closed-loop  systems.  But  first  it  has  to  be  clarified  what 
is  meant  by  det[X  +  'yTQ{s)]  for  the  infinite-dimensional  operator  I  +  ')TQ{s). 

From  Assumption  (*)  we  know  that  ^(s)  G  3§i{X)  for  some  s  G  p{A).  Then  it  is  simple  to  show  that 
0(s)  G  ^i(X)  for  every  s  G  p(A)  [9].  Also  T  G  SS(fC)  implies  TQ{s)  G  ^i(X).  One  can  now  define  [10]  [9] 

dei[I  +  yTQ{s)]  :=  (l  +  7An(s)). 


where  {^n{s)} are  the  eigenvalues  of  G{s). 

On  the  other  hand,  the  boundedness  of  the  operators  B,  C,  T  together  with  Assumption  (*)  imply  that  (a)  A 
and  71“^*  =  A  —  ^BTC  are  defined  on  the  same  dense  domain  '3,  (b)  p{A)  H  is  not  empty,  (c)  for  all 

s  G  p{A)  we  have  yBTC{sI  —  A)~^  G  This  allows  us  to  introduce  the  perturbation  determinant  [11] 

A_4d/_4(s)  :=  dei[{sl  —  A'^’^){sl  —  A)~^] 

=  dei[l  +  yBTC{sl  -  Ay^]  =  det[l  +  yTg{s)] 


which  is  analytic  in  p{A)  H  ^(-4“^*)  [see  Lemma  Al].  In  fact  A^d /_4(s)  is  the  equivalent  of  the  fraction  in  (1)  for 
systems  with  unbounded  infinitesimal  generators.  We  are  now  ready  to  state  an  extended  form  of  the  argument 
principle  for  such  systems.  The  following  theorem  makes  use  of  the  formula  [11] 

^  In  A_4d/_4(s)  =  tr[(sX  —  —  (sX  —  ^)“^]  for  all  s  &  p{A)  p{A'^')  (4) 

to  relate  det[X  +  7-^^('5)]  to  the  eigenvalues  of  A  and  A'^^  inside  D. 

Theorem  1:  If  det[X  +  'yTQ{s)\  0  for  all  s  G  D, 


C^OidetiX-f  7X'e;(s)] 


sG'D 


=  tr 


[(si  -  Ay^ds 


—  tr| 


- —  /  (si  —  A)  ^ds 
2ttj' 


t'S 


=  —  (number  of  eigenvalues  of  A‘^^  in  C“'') 

-f  (number  of  eigenvalues  of  A  in  C“''), 


where  D  is  the  Nyquist  path  shown  in  Figure  2  that  does  not  pass  through  any  eigenvalues  of  A. 

Proof:  See  Appendix.  ■ 

Remark  1:  Theorem  1  relies  on  the  fact  that  under  Assumption  (*)  both  {si  —  A)~^  and  {si  —  are 

compact  operators,  which  means  the  infinitesimal  generators  A  and  A'^^  have  discrete  spectrum  (i.e.,  their  spectrum 
consists  entirely  of  isolated  eigenvalues  with  finite  multiplicity).  Then  V  =  —  A)~^ds  is  the  group- 

projection  [12]  [9]  corresponding  to  the  eigenvalues  of  A  inside  D,  and  tr[X]  gives  the  total  number  of  such 
eigenvalues  [13].  Similarly  tr[— 2:^  ~  gives  the  total  number  of  eigenvalues  of  71“^’  in  J).  ■ 


APRIL  17,  2006 


4 


As  a  direct  consequence  of  Theorem  1  we  have  the  following. 

Theorem  2:  Assume  p+  denotes  the  number  of  eigenvalues  of  A  inside  C’*'.  For  S  taken  as  above,  the  closed- 
loop  system  is  stable  iff 

(a)  det[l  + -iTg{s)]  ^  0,  V  s  G  S), 
and 

(b)  c(0;det[X  +  7.Fg(s)]|^^^)  =  p+. 


B.  The  Eigenloci  Method 

The  setback  with  the  method  described  in  the  previous  paragraph  is  that  to  show  C  C“,  = 

A  —  B"f!FC,  for  different  values  of  7,  one  has  to  plot  det[X  +  7-^^(s)]  for  each  7.  Note  that  this  includes 
having  to  calculate  the  determinant  of  an  infinite  dimensional  matrix.  This  motivates  the  following  eigenloci 
approach  to  Nyquist  stability  analysis,  which  is  very  similar  to  that  performed  in  [4]  for  the  case  of  time-periodic 
systems. 

Let  {^n(s)}„g2  constitute  the  eigenvalues  of  TQ{s).  Then 

Z  det[X  +  7.F0(s)]  =  Z  J]  (1  +  7A„(s)).  (5) 

Recall  that  TQ{s)  G  S§i{X)  for  every  s  G  p{A).  This,  in  particular,  means  that  TQ{s)  is  a  compact  operator  and 
thus  its  eigenvalues  A„(s)  accumulate  at  the  origin  as  |n|  — >  00  [14].  As  a  matter  of  fact  one  can  make  a  much 
stronger  statement. 

Lemma  3:  The  eigenvalues  A„(s),  s  G  S,  converge  to  the  origin  uniformly  on  S. 

Proof:  See  Appendix.  ■ 

Take  the  positive  integer  to  be  such  that  |A„(s)|  <  e,  s  G  S,  for  all  \n\  >  N^.  Let  us  rewrite  (5)  as 

Z  det[I  +  -iTQ{s)]  =  Z  []  (i  +  7A„(s))  +  Z  J]  (i  +  7A„(s)) 

\n\<N^  |n|>Afe 

=  ^  ^  (l+7'^n('S))  +  ^  ^  (1  +  7'^n('S))-  (6) 

InlGAf^  |n|>tV5 

It  is  clear  that  if  I7I  <  f  then  for  |n|  >  we  have  |7A„(s)|  <  1,  and  l  +  7A„(s)  can  never  circle  the  origin  as  s 
travels  around  S.  Thus  for  I7I  <  ^  the  final  sum  in  (6)  will  not  contribute  to  the  encirclements  of  the  origin,  and 
hence  we  lose  nothing  by  considering  only  the  first  N^;  eigenvalues.  There  still  remain  some  minor  technicalities. 

First,  let  denote  the  disk  |s|  <  e  in  the  complex  plane.  Then  said  truncation  may  result  in  some  eigenloci 
(parts  of  which  reside  inside  D^)  not  forming  closed  loops.  But  notice  that  these  can  be  arbitrarily  closed  inside 
D^,  as  this  does  not  affect  the  encirclements  [4]. 

The  second  issue  is  that  for  some  values  of  s  G  D,  d^G{s)  may  have  multiple  eigenvalues,  and  hence  there 
is  ambiguity  in  how  the  eigenloci  of  the  Nyquist  diagram  should  be  indexed.  But  this  poses  no  problem  as  far 
as  counting  the  encirclements  is  concerned,  and  it  is  always  possible  to  find  such  an  indexing;  for  a  detailed 
treatment  see  [3]. 

Let  us  denote  by  indexed  eigenloci  that  make  up  the  generalized  Nyquist  diagram.  [To  avoid 

confusion  we  stress  the  notation;  A„(s)  is  the  n*  eigenvalue  of  TQ{s)  for  a  given  point  s  G  D,  whereas  A„  is 
the  n*  eigenlocus  traced  out  by  A„(s)  as  s  travels  once  around  D.]  From  (6)  and  the  above  discussion  it  follows 
that 

C'(0;det[XH-7.F5(s)]|^^^)  =  ^  c(--;A„) 

which  together  with  Theorem  2  gives  the  following. 

Theorem  4:  Assume  p+  denotes  the  number  of  eigenvalues  of  A  inside  C+.  For  S  and  as  defined  previously, 
the  closed-loop  system  is  stable  for  I7I  <  ^  iff 
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(a) 
and 

(b) 


X]|n|<7Ve  ~  =  P+- 


V.  An  Illustrative  Example 

Consider  the  system  defined  on  the  interval  a;  €  [0, 27r]  and  governed  by  the  Partial  Differential  Equation  (PDE) 
dtip{t,x)  =  d‘^tp{t,x)  —  ^cos{x)  ip{t,x)  +  'ilj{t,x), 
with  7  G  C  and  periodic  boundary  conditions 

0)  =  'tp{t,2TT),  dx'ip{t,0)  =  dxi’{t,2'K). 

Let  us  rewrite  this  system  in  the  form  of  a  PDE  described  by 

dt^{t,x)  =  d‘l'ip{t,x)  +  ^{t,x)  +  u{t,x), 

y{t,x)  =  tp{t,x),  (7) 

placed  in  feedback  with  the  function 

'yF{x)  =  7Cos(a;). 

The  problem  is  now  in  the  general  from  discussed  in  Section  II  and  can  be  considered  as  a  differential  equation 
on  X  =  L^[0, 27r];  A  =  +  1  and  is  defined  on  the  dense  domain 

^  =  {4>  £  L^[0,  27r]  I  (j),  ^  absolutely  continuous,  G  L^[0, 27r],  (/)(0)  =  (/)(27r),  ^(0)  =  ^(27r)}, 

ax  dx^  dx  dx 

B  and  C  are  the  identity  operator,  and  F  =  cos(a;). 

We  take  an  extra  step  and  use  a  similarity  transformation  to  put  the  problem  in  an  equivalent  form  that  is 
more  familiar  to  us  from  multivariable  linear  systems  theory.  Let  ^  be  the  transformation  that  takes  the  function 
(j){x)  G  E^[0,  27r],  (j){x)  =  X]n6Z  4>n  6-^"®,  to  its  Eourier  series  coefficients  col[-  •  •  ,  (/)o,  (/>i,  •  •  •  ]  G  Then  it 

is  simple  to  show  that  A,  B,  C,  and  F  have  the  following  (bi-infinite)  matrix  representations 


,  B  =  C  = 

F  = 


0  1 
1  0 


Since  S(g1)  =  {— n^  +  l,n  G  Z},  then  Assumption  (**)  is  satisfied.  Eor  any  s  ^  E(A)  we  have  (s  —  ^)  ^  = 
diag{---  ,  5+7231,- ••}■  Thus  E„gN =  E„ezL+J^-il  <  o®-  Hence  G 

and  Assumption  (*)  is  satisfied. 

Notice  that  the  open-loop  system  is  unstable.  Next  we  demonstrate  that  by  plotting  the  eigenloci  one  can  read 
off  from  this  plot  the  stability  of  the  closed-loop  system  for  any  value  of  7  G  C. 

A  =  0,  0, 1  are  the  eigenvalues  of  A  inside  2),  hence  p+  =  3,  and  we  need  three  counter-clockwise  encirclements 
of  — 1/7  to  achieve  closed-loop  stability.  As  can  be  seen  in  Eigure  3(b)  and  its  blown-up  version  (c),  one  possible 
choice  would  be  to  take  — 1/7  to  be  purely  imaginary  and  —  0.2j  <  — 1/7  <  0.2 j.  Clearly  such  — 1/7  is  encircled 
three  times  by  the  eigenloci. 


VI.  Conclusions 

We  develop  an  extension  of  the  Argument  Principle  and  the  Nyquist  Stability  Criterion  that  is  applicable  to 
systems  with  infinitesimal  generators  that  are  unbounded  operators  with  discrete  spectrum  and  whose  resolvent 
operator  is  trace-class.  This  theory  can  be  used  to  verify  the  stability  of  spatially  extended  systems  and  those 
governed  by  partial  differential  equations,  as  demonstrated  in  an  example. 
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Fig.  3.  The  Nyquist  path  X>;  The  Nyquist  plot;  Blown-up  version  of  the  center  part  of  the  Nyquist  plot. 


VII.  Appendix 

To  prove  Theorem  1  we  need  the  following  lemma. 

Lemma  Al:  For  s  S  p{A),  det[X  +  ^!FQ{sy\  is  analytic  in  both  7  and  s. 

Proof:  For  s  €  p{A),  'yJ-Q{s)  €  ^i(X).  Also  —  A)~^B  is  clearly  analytic  in  both  7 

and  s  for  s  G  p(A).  Then  it  follows  from  [9,  pl63]  that  det[I  +  too  is  analytic  in  both  7  and  s  for 

s  e  p(A).  m 


Proof  of  Theorem  7:  Consider  any  point  s  in  D.  Since  J)  does  not  pass  through  any  eigenvalues  of  .4,  s  G  p{A) 
and  thus  ^TC{sI  —  A)~^B  G  P§i{X)  by  Assumption  (*).  Then  from  [10],  (l  +  ^TC{sI  —  A)~^B)  ^  exists 
and  belongs  to  SSifX')  iff  det[X  +  ^TC{sX  —  A)~^B]  f  0,  which  is  satisfied  by  assumption.  Applying  an 
operator  version  of  the  matrix  inversion  lemma  to  (X  +  ^TC{sI  —  A)~^B')  ,  we  conclude  that  {sX  —  = 

{sX  —  A  +  B^TC)~^  G  PS{X)  and  thus  s  G  Therefore  D  is  contained  inside  p{A)  H  p{A'^^). 

Let  the  path  €  be  that  traversed  by  det[X+7iF^(s)]  as  s  travels  once  around  D.  By  Lemma  Al,  det[X+7iF^(s)] 
is  analytic  in  s,  and  if  det[X  +  ^TQ{s)]  0  on  £)  we  have 

C(0;det|l  +  7«(.)l|.,^)  =  XJL 

^  f  g;det[X+7.Xg(s)] 
det[X  +  7X'^(s)] 

1  f  d 

=  — —  /  tv[(sX  —  —  (sX  —  A)^^]ds,  (Al) 

27rj  7s, 

where  we  have  used  (4)  in  the  last  equality.  Notice  that  because 

{sX  -  A^^Y^  -  (sX  -  A)-^  =  -7(sJ-^)-iSX-C(sX-X')-^  G  ^i(X)  (A2) 

for  all  s  G  p(4l)np(^‘^*),  the  right-hand  side  of  (Al)  makes  sense  and  is  finite.  (A2)  also  gives  that  {sX—A'^^)~^  G 
,!^l(X). 

On  the  other  hand,  since  {sX  —  A)~^  and  {sX  —  both  belong  to  7^i(X)  C  i^oo(X),  their  spectra  consist 

entirely  of  isolated  eigenvalues  with  no  finite  accumulation  point  [13,  pi 87].  By  Assumption  (**)  the  real  part 
of  the  essential  spectrum  of  A  is  —00.  Now  from  (A2)  we  have  that  {sX  —  —  {sX  —  A)~^  G  ,^00 (X)  and 

thus  A^^  has  the  same  essential  spectrum  as  A  [13,  p244].  Therefore  the  path  S)  encloses  a  finite  number  of  the 
eigenvalues  of  A  and  Thus  in 

[(sX-A'^Y^ds  -  [(sX-A)~^ds 

27rj  2ti] 
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each  term  is  a  finite-dimensional  projection  [9,  pll,  pl5].  Taking  the  trace,  from  [13]  it  follows  that 


tr 


[ {si  —  ^ds  —  tr  [ {si  —  A)  ^ds 


(A3) 


is  equal  to  the  number  of  eigenvalues  of  ^  in  2)  minus  the  number  of  eigenvalues  of  in  2),  where  2)  is  the 
(clockwise)  Nyquist  path  and  is  taken  arbitrarily  large  to  enclose  C+.  Finally  (A3)  and  (Al)  together  give  the 
required  result.  ■ 


Proof  of  Lemma  3:  For  s  G  2)  C  p{A),  dei\I  + ')TQ{s)]  is  analytic  in  both  7  and  s  by  Lemma  Al.  The  proof 
now  proceeds  exactly  as  in  [4,  pl40]  and  is  omitted.  ■ 
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